Electron-hole liquid in the hexaborides 
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We investigate the energetics of the electron-hole liquid in stoichiometric divalent metal hexa- 
borides. The ground state energy of an electron-hole plasma is calculated using RPA and Hubbard 
schemes and compared to the binding energy of a single exciton. Intervalley scattering processes 
play an important role in increasing this binding energy and stabilizing a dilute Bose gas of excitons. 
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The remarkable discovery of the high-temperature 
weak ferromagnetism in La-doped CaBg, SrBg, and BaBg 
has opened a new page in the physics of magnetism pi . 
In our previous work J2| we attributed this effect to an 
unusual ground state of undoped divalent hexaborides. 
This so-called excitonic insulator is characterized by a 
condensation of bound electron-hole pairs (excitons). An 
excitonic instability in narrow gap semiconductors or 
semimetals was predicted Q and studied theoretically in 
the mid-sixties [||. However, its occurrence in any real 
compound is still controversial. Band structure calcula- 
tions [gU^l predict a small direct overlap in divalent-metal 
hexaborides between a boron-derived valence band and 
a cation-derived conduction band at three equivalent X- 
points in the cubic Brillouin zone. This feature, together 
with absence of direct electric-dipole transitions between 
the two bands, is extremely favorable for electron- hole 
pairing and leads to the excitonic instability. Weak fer- 
romagnetism develops, then, in a triplet excitonic insula- 
tor due to spontaneous time-reversal symmetry breaking 
under doping [^|j. 

There are many open questions in the physics of the 
excitonic ferromagnetism. Some of them were raised re- 
cently in Refs. II^,^]. In the present work we, however, 
want to shift emphasis from the (un)doped excitonic state 
to general properties of electron-hole {e-h) liquids in the 
hexaborides. Charge conservation does not fix the num- 
ber of e-h pairs, which in thermal equilibrium depends on 
the e-h chemical potential (band overlap). In the original 
theory of the excitonic insulator the band gap (overlap) 
Eg > {Eg < 0) was considered to be a free parame- 
ter that changes continuously through the value Eg = 0. 
Subsequently, investigation of optically pumped electrons 
and holes in semiconductors showed |lO|-[l4| that a first 
order transition between two states, one with a substan- 
tial band gap and the other with a substantial band over- 
lap, will occur and that smaller values of \Eg\ lie in the 
unphysical intermediate region. The details of such tran- 
sition are very sensitive to the actual band degeneracies 
and anisotropics. Here, we study possible scenarios for 
the transition from a semiconducting state into a metal- 
lic e-h liquid in CaBg, including the appearance of a free 
exciton gas. We find that a novel mechanism of interval- 



electron mass [g: 
tion band) and ml^ 



ley scattering of excitons is important in stabilizing an 
intermediate gas phase. 

Electrons and holes in CaBg have the following val- 
ues of the effective masses measured in units of the bare 
= 0.504, m^ = 0.212 (conduc- 
2.17, m^ = 0.206 (valence band). 
These values agree with the early results using the muffin- 
tin approximation Iq]. First, we retain only the dom- 
inant intravalley scattering processes shown in Fig. ^a, 
when an electron (hole) scatters between states near a 
single AT-point. In the small-g limit the scattering ma- 
trix element is given by a screened Coulomb potential 
Vq = Ane'^/Kq'^ {k is a static dielectric constant). We 
define an e-h chemical potential as a sum of two in- 
dividual potentials /i — fie -{- fJ-h- The relation to the 
band model is made hy fj, — —Eg, i.e. the band gap 
is equivalent to the chemical potential. Natural units 
for energies and lengths are effective Rydberg E^ and 
Bohr radius Qx- Ex = TO*e*/2K^ = e^ jl^ax, where re- 
duced mass TO* = moeTnohl {tnoe + TiT'oh) is determined 
by optical masses: 3/too = Ijm^ -f 1/to". The e-h pair 
density n is characterized by a dimensionless parameter, 
r, = (3/47rna3)i/3. 
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FIG. 1. Electron-hole scattering processes in the hexa- 
borides. Operators a and 6 correspond to conduction and 
valence bands, respectively. Figure (a) shows a dominant in- 
travalley vertex. Figures (b) and (c) correspond to intervalley 
scattering. 

At high densities, small r^, we use the random phase 
approximation (RPA) |ll|. Strictly speaking, a dense 
metallic e-h liquid can transform into an excitonic in- 
sulator [pP]. This assumption was a key starting point 
in our explanation of the high-temperature weak ferro- 
magnetism in the hexaborides B. However, the ground 



state energy correction from the excitonic instability in 
a semimetal is of order i^ jsp and is small (A <C Ef), so 
we neglect it. The kinetic energy per e-h pair is given by 
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where a — (4/97r)^/"^ and v (= 3) is the number of valleys. 
The exchange energy for anisotropic bands is 



E, 



cxch 



2-KarsV^I'^ 



(/)(to^/to||) + 0(m;|:/m|[) 



(x) — X 



]^/g arcsin-v/T" 



VT^^ 



(2) 



The correlation energy describes the remaining contribu- 
tion to Err s 
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where n*(g,w) is the irreducible polarization operator 
and U"{q,uj) = X^i" n°(g,u;) is a sum of (anisotropic) 
RPA-polarizabilities for each species of electrons or holes. 
Substitution U*{q,Lu) = n°(g,w) in Eq. (|) gives the 
RPA-expression for Ec. An approximate way of treat- 
ing the higher order exchange corrections was considered 
by Hubbard [|3|. His expression generalized to the mul- 
ticomponent plasma is 
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with f{q) — 0.5q^/{q'^ + kp). We also change an uj- 
integration from real to imaginary axis, which avoids 
a difficulty related to a plasmon pole in U'^{q,uj). Nu- 
merical results for the ground state energy are shown in 
Fig. y. The band degeneracy improves substantially the 
accuracy of the RPA, because corrections to the RPA di- 
agrams have extra smallness in l/2z/. The minimum of 
the ground state energy is reached at r^ — 0.92 with a 
minimum value E^l"^ = -1.55E.^{-l.51E^) in the RPA 
(Hubbard) scheme. The use of RPA is justified by a small 
value of Ts at the minimum, which corresponds to a dense 
plasma with strong screening and small corrections from 
multiple e-h scattering. Such corrections become signifi- 
cant at Ts ^ 3. The ground state energy is also known in 
the limit rs — > oo, where it approaches the binding en- 
ergy of a single exciton E^- (Here, we disregard possible 
formation of exciton molecules.) 

The presence of a local minimum at metallic densi- 
ties for Egs(rg) has an important effect on transforma- 
tion from a semiconductor to a semimetal |10|-p2[. The 
two possible scenarios are shown schematically in Fig. H. 
In the first case, curve (a), the local minimum is also 
the absolute one: -Bg.s.(''sA) < E^- The pair chemi- 



M 



Err 



dE, 



dn 



At the extremal point the second 
term is zero and /i = E'^l^. Therefore, when the band gap 
decreases to Ec — |£^™™|, a first-order metal-insulator 
transition takes place. The number of e-h pairs jumps 
from zero to n{rsA)- Smaller densities with r^ > Tsa 
correspond to unstable states. In optically pumped semi- 
conductors, where a number of carriers is fixed instead of 
a chemical potential, this effect leads to e-h droplet con- 
densation ||l^]. In the second case, curve (b), the exciton 
energy lies below the metallic minimum. Semiconduct- 
ing state becomes unstable at Eg = \Ee\ and transforms 
into a low-density exciton gas. If Eq is further reduced, 
a first order transition of gas-liquid type takes place be- 
tween two states Bi and B2, which have same pressure 

,2 9Bg.. 



P ~ n^ "an" • ^^^ *"^^^ (^) ^^ believed to be realized 
in isotropic one-component e-h plasma, whereas in many 
semiconductors band degeneracies and anisotropics favor 
the case (a) bypassing a free exciton gas state ||ll| , [l^ . 
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FIG. 2. The ground state energy per pair in units of E^c 
for e-h liquid in CaBe. Solid line is the RPA and dashed line 
is the Hubbard result. Dot-dashed line is the RPA including 
correction from the intervalley scattering. Thin horizontal 
lines denote position of the Ai exciton level in zeroth (upper) 
and first (lower) order approximations in the intervalley scat- 
tering. Self-consistent mechanism proposed in the text pushes 
the exciton energy further down below the local minimum. 

To check which of the two scenarios occurs in the hex- 
aborides we now compare the above value of i?""" to the 
binding energy of a single exciton. The simplest esti- 
mate for the exciton energy is —Ex- It is obtained by 
substituting an isotropic Is-type wave function into the 
Schrodinger equation. This estimate predicts the same 
binding energy for excitons formed by an electron and 
a hole from same or from different valleys. We improve 
this result by using a simple variational ansatz appropri- 
ate to the cylindrical symmetry of the bands near each of 
the AT-points similar to the treatment of shallow impu- 
rity states [|Ij| . The binding energy of an electron and a 
hole from the same valley is E^, = —1.12Ex {{a/xbx) 7^ 0), 



cal potential is related to the ground state energy by: whereas an electron and a hole from different valleys have 



E'^ = -IME^ ((flx^x') "^ 0)- The former case with a 
more anisotropic reduced mass tensor is favored since 
lower dimensionahty gives, as usual, a stronger binding. 
Thus, band anisotropy favors a particular type of exci- 
tons diagonal in the valley index. The chemical potential 
of the metallic e-h liquid is below the exciton binding en- 
ergy in this approximation. Hence, a dilute exciton gas 
is unstable and scenario (a) is predicted instead. 
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FIG. 3. Schematic form of the density dependence of 
the ground-state energy per pair in the two scenarios for 
metal-insulator transition in a semiconductor. 

The dominant term approximation (Fig. |l|a) used 
above leaves threefold degeneracy of the excitonic level 
corresponding to e-h pairing at three X-points. Coher- 
ence between different excitons is established by inter- 
valley scattering processes shown in Fig. Mj. Instead 
of a single hydrogenic equation, one must solve now a 
system of three coupled integral equations for a three- 
component excitonic wave function jpi, (r) . The scattering 
matrix element Qq includes a large momentum transfer 
on Q = &/v2, where b is an elementary reciprocal lat- 
tice vector. Consequently, we can neglect its momentum 
dependence and estimate 



E 



47re^ 



IQ + GP 



F^°(X,X')F°'g(^'.^ 



(5) 



Here, summation is over reciprocal lattice, Fq''{X,X') 
is a form factor 0] of two conduction band states at 
points X and X'. The absence of the dielectric con- 
stant in the Fourier transform of the Coulomb potential 
reflects a lack of lattice screening for these processes. An- 
other vertex, which has the same strength gq but does 
not contribute to the electron-hole pairing, is shown in 
Fig. nr. Several other intervalley vertices correspond to 
processes, when an electron transforms into a hole. How- 
ever, they correspond to higher momentum transfers and 
we neglect them. Generally, the form factors in Eq. (0) 
reduce the matrix element compared to the amplitude 



of a Coulomb potential by a factor of ~ 2-5 and can 
also change the sign of gq. Nevertheless, in the absence 
of necessary numerical band structure results we use a 
somewhat optimistic estimate \g\ « Aire'^ / Q'^ . Returning 
back to a spatial form of the Schrodinger equation we 
obtain a contact-like interaction of excitons in different 
valleys with a strength g. 

Since the total momentum of excitons vanishes, eigen- 
statcs of the Schrodinger equation are classified accord- 
ing to the irreducible representations of the cubic point 
group Oh- Treating the intervalley scattering pertur- 
batively, we find E^i = Ee — 2X for a nondegenerate 
exciton state with Ai symmetry ^^ ~ '0B(r)(l, 1, 1) 
and Ee2 = E^ + X for a degenerate doublet with E 
symmetry: ipi^^ - Vb(^)(1, e2"/^ e'*''*/^) and tpP ~ 
i/'s(?')(l,e~^'^*/'^,e~'*'^*/'^), where Vs('') is a Is hydro- 
genic wave function. The coupling constant is A = 
g\tpB{0)\'^ and, therefore, choice of the lowest energy ex- 
citon depends on the sign of g. If interaction between 
excitons is attractive {g > 0) the symmetric ^i-state 
has lower energy. If it is repulsive, then the i?-doublet 
is more stable. One can easily estimate A using unper- 
turbed ipBif)- A = 0.15Ex. We have also calculated the 
effect of intervalley vertices Figs. 0b and np on the ground 
state energy of the metallic e-h liquid shown by dashed 
line in Fig. S. This correction is much smaller than a 
change of the exciton energy and does not exceed 3%. 
Qualitatively, such a difference is explained by different 
orders of the two corrections: for the metallic plasma it 
is a second-order effect, while the shift of exciton energy 
is a first order effect. Another effect of lifting threefold 
degeneracy of excitons in different valleys is suppression 
of multi-exciton molecules. 

As one can see from Fig. y the exciton energy is still 
above the ground state energy of the metallic phase, 
though the two energies move closer to each other. How- 
ever, it appears that our estimate for the intervalley scat- 
tering effect may be too conservative. The origin of extra 
enhancement is similar to a mechanism of anomalously 
large hypcrfinc splitting of donor states in Si proposed 
by Kohn and Luttingcr ||l^. The actual exciton energy 
is shifted below prediction of the effective mass theory 
because of both the intervalley scattering effect and a so 
called central cell correction [^. Contrary to a naive 
point of view, related changes in the amplitude ■0(0) 
are determined by a long distance Schrodinger equation 
rather than by an exact short distance e-h Hamiltonian. 
Since the actual Ee is not an eigenvalue of the effective 
mass equation, there is no solution for this energy which 
satisfies both boundary conditions at r = 0,cxd. There- 
fore, actual exponentially decaying wave-function devel- 
ops a singularity at short distances. We find that for a 
moderate 20% shift of E^ the probability to find electron 
near hole |V'(0)P is enhanced by a factor of 4 compared 
to |V'b(0)P. This effect significantly increases A, espe- 



daily for the Ai-singlet [Q. As a result, the excitonic 
level in the hexaborides must lie well below the metallic 
minimum of -Eg.s.(^s)- 

We have considered stability of different phases of the 
e-h liquid in the hexaborides. They include a semicon- 
ducting state with no carriers, a dilute Bose gas of ex- 
citons and a dense electron-hole liquid. The latter can 
be clearly distinguished from the other two states in in- 
frared optical conductivity by its large Drude peak. If an 
excitonic instability develops in a dense e-h liquid, (j{uj) 
must also have an edge-type singularity at the excitonic 
gap. Semiconducting and free exciton gas states, on the 
other hand, have no significant features in (t{uj) because 
of the absence of direct optical transitions between the 
two bands. Which of the three states really occurs de- 
pends on the band gap parameter Eq- Note, that in the 
scenario (a), see Fig. |[ a first order transition between 
a semiconducting state and a dense plasma takes place 
at a positive Eq, which always corresponds to a semi- 
conductor in a single electron picture. Small values of 
\Eg\ ^ Ex = 0.08 eV found in the band structure cal- 
culations [^J|] indicate that the hexaborides can be close 
to an instability of semiconducting state, which, as we 
argued for CaBg, transforms into a dilute exciton gas. 
The other possible instability in the hexaborides, if Eq 
is varied on experiment, is a gas-liquid transition between 
a dilute Bose gas of excitons and a dense e-h plasma. 

We suggest that some of the above phase transforma- 
tions could be induced in the hexaborides by applying 
hydrostatic pressure, which changes the band gap Eq. 
Applied uniaxial stress can further lift the degeneracy 
between different valleys and, thus, reduce significantly 
an energy of the metallic phase [[ll|jl^ ]. Our main result, 
stability of a dilute exciton gas in CaBg either at nor- 
mal conditions or under pressure, provides a new way to 
look at the Bose condensation of excitons. Another in- 
trinsic mechanism for gap variations can also come from 
impurity doping, in particular from La- substitution for 
a divalent element. Therefore, doped hexaborides can 
differ from undoped compounds not only because of un- 
equal number of electrons and holes, but also in terms of 
an appropriate starting picture: a Bose gas of excitons 
or a dense e-h liquid. In the latter case screening effects 
must be quite significant due to a high density restriction 
for e-h plasma: rg ^ 1. 
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Central cell corrections to the exciton energy come in part 

because of a lack of lattice screening of e-h interaction 

at distances smaller than lattice constant. We estimate 

this effect as lS.Ee ~ ~Q.2Ex for the Ai-state. Note that 

similar shift is also present for the ground state energy 

of e-h plasma with r^ ~ 1. Hence, central cell correction 

alone cannot explain stability of a dilute exciton gas. The 

i5-states have vanishing probability to find an electron at 
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them. 
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